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ABSTRACT 

Event shapes have long been used to extract information about hadronic 
final states and the properties of QCD, such as particle spin and the run- 
ning coupling. Recently, a family of event shapes, the angularities, has 
been introduced that depends on a continuous parameter. This additional 
parameter-dependence further extends the versatility of event shapes. It 
provides a handle on nonperturbative power corrections, on non-global 
logarithms, and on the flow of color in the final state. 
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1 Introduction 



Event shapes are generalizations of jet cross sections that describe the distribution 
of radiation in the final state. Because no individual hadrons are observed, event 
shapes are infrared-safe. This renders the main features of event shapes, which are 
determined by the underlying hard scattering, perturbatively calculable (for recent 
summaries see Refs. [1, 2]). Event shape observables are therefore an ideal testing 
ground of QCD. For example, they allow the determination of the gluon spin [3], of 
QCD color factors [4], and precise measurements of the running coupling [5]. 

Nevertheless, event shape observables retain sensitivity to long-distance hadroniza- 
tion effects. For average values of e + e~ event shapes, these nonperturbative cor- 
rections manifest themselves as more-or-less additive corrections, typically propor- 
tional to the inverse power of the center-of-mass (cm.) energy, 1/Q. Differential 
distributions, which give important additional information about the dynamics of 
QCD, however, require a more involved convolution with functions that parameterize 
hadronization effects (for reviews see for example Refs. [1, 6]). 

Below we will illustrate how the introduction of a parameter-dependence into 
an event shape observable can considerably extend the information that can be ex- 
tracted from its measurement. We discuss a set of event shapes, the so-called class of 
angularities [7, 8, 9], that depends on a continuous real parameter. This parameter- 
dependence allows to test nonperturbative features of event shape distributions [10], 
to control large non-global logarithmic corrections [7, 8], and, in the case of hadronic 
collisions, to obtain information about the underlying color flow [11]. 

We begin with a necessarily incomplete review of a prominent event shape, the 
thrust [12, 13], and introduce the family of angularities. We then derive the scaling 
rule for nonperturbative power corrections to distributions of angularities. Next, 
we illustrate how the correlation of angularities with non-global observables, that is, 
with observables that are designed to measure radiation into only part of phase space, 
can be used to dampen the sensitivity of the perturbative calculation to soft-gluon 
radiation at wide angles away from the point of interest. In Section 5 we propose an 
extension of the definition of angularities to hadronic collisions, and show how such an 
observable may be used to extract information about the distribution of color in the 
underlying hard scattering. We end with a summary and point out future directions. 



1 



2 The Thrust and Its Extension to the Family of 
Angularities 



2.1 Determination of Gluon Spin via the Thrust 

One of the prime examples of an event shape observable in e + e~ collisions is the 
thrust [12, 13] T: 

T = ———maxJ2\Pi-n\ , (1) 

where the sums are over all particles in the final state N, and n is a unit vector 
whose direction is called the thrust axis when T is maximal. For massless partons, 
this definition can be written equivalently as 

r = l-r=i£ p ± ,ie-M, (2) 
^ ieN 

where p±^ and rji denote the transverse momentum and pseudorapidity rji = In cot(0j/2), 
respectively, of particle i, with respect to the thrust axis. The maximization is implied. 
Here and below, we treat all partons as massless, parton masses induce corrections 
which have been examined in Ref. [14]. 

The thrust measures how pencil-like a particular event is. In the limit of two 
back-to-back jets, the value of T is 1, whereas its minimum value of 1/2 corresponds 
to a completely homogeneous spherical event. Radiative corrections shift the value 
of the lowest order cross section away from T — 1. The distribution of radiation, and 
thus the differential cross section da/dT, depends on the spin of the radiated partons. 
The first-order QCD cross section is given by 



1 da(qqg) 2a s 



dT 3tt 



2(3T 2 -3T + 2) /2T- 1\ 3(3T-2)(2-T) 
T(l-T) 11 V 1 - T ) 1 — T 



(3) 



whereas the corresponding cross section where a scalar "gluon" is radiated gives a 
distinctly different distribution, 



1 da(qqS) a & 



dT 3tt 



2 In ( 2T - l \ | (4~3T)(3T-2) 



V 1 -T J 1 — T 



(4) 



This sensitivity to the spin of the radiated partons was used as one of the first tests 
of QCD [3]. The definitions of a multitude of other e + e~ event shapes measuring 
various aspects of angular distributions in the final state can be found, for example, 
in Ref. [1]. 
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2.2 The Family of Angularities 



We now extend the definition of the thrust, Eq. (4), by introducing a dependence on 
a continuous, real parameter a, which allows us to study a whole class of event shapes 
simultaneously. We define the family of angularities by weighing the final state with 
the function [7, 8, 9] 

Ta(N) = ± £ P±,i e ~ Mil ~ a) > ( 5 ) 

°* ieN 

where the direction of particle i is again measured relative to the thrust axis. The 
parameter a varies between a = 2 (the weight function (5) then selects two infinitely 
narrow, back-to-back jets) and a = — oo (the fully inclusive limit). The value a = 
corresponds to 1 — T, with T the thrust (1), and a = 1 corresponds to the jet 
broadening [15]. 

In the two-jet limit, r a approaches 0, and the differential cross section da/dr a 
receives large corrections in ln(l/r a ) due to soft gluon radiation which have to be 
resummed for a reliable theoretical prediction. This resummation has been performed 
to all logarithmic orders at leading power for a < 1 in Ref. [8]. For a ~ 1, recoil 
effects have to be taken into account, as was pointed out for the broadening (a = 1) 
in Ref. [16]. All equations below are therefore valid in the range a < 1, where the 
value of a is not too close to I. 

We quote the result of the resummation of large logarithms of r a in Laplace mo- 
ment space: 

<r{i/,Q,a)= dr a e . (6) 

Jo dr a 

Logarithms of l/r are transformed into logarithms of v. At next-to-leading loga- 
rithmic (NLL) order, we find 



— o (v, Q,a) = exp j 2 / *t / (a s (p ± )) _ i) 



We see that for small r a , the cross section factorizes into contributions from two 
independently radiating jets, resulting in the Sudakov exponent above, where the 
factor of two stems from the two equal jet contributions. Here, A(a s ) and B(a s ) 
are well-known anomalous dimensions that are independent of a. They have finite 
expansions in the running coupling, A(a s ) = Xi™=i (a s /7r)", and similarly for 
B(a 8 ), with the coefficients A« = C F , = -3/2 C F , A^ = 1/2 C F [C A (67/18 - 
7r 2 /6) — 10/9 T F Nf], at NLL accuracy, with the color factors Cp = 4/3, Ca = 3, 
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T F = 1/2, and Nf denotes the number of flavors. At a = we reproduce the NLL 
resummed thrust cross section [17]. The explicit NLL expression in transform space 
can be found in Ref. [10], and a formula for the resummed cross section valid to all 
logarithmic orders at leading power is given in Ref. [8]. 



3 Universality of Power Corrections 

As already mentioned in the introduction, the resummed expression (7) is plagued by 
sensitivity to long-distance effects. This sensitivity manifests itself as an ambiguity 
in how to deal with the not well-defined integral over the running coupling, and the 
validity of the perturbative approach breaks down at r a ~ A QCD /(5. There exist a 
variety of approaches (see for example Ref. [6]), but in all cases one has to supplement 
the perturbative calculation with additional information. In the present case, as we 
will see below, the knowledge of nonperturbative (NP) corrections for only one specific 
value of a suffices to determine unambiguously the full distribution, including power 
corrections, for all other values of a (a < 1). 

Due to the quantum mechanical incoherence of short- and long-distance effects, 
we can separate the perturbative part from the NP contribution in a well-defined, 
although prescription-dependent manner. Following Refs. [18, 19], we deduce the 
structure of the NP corrections by a direct expansion of the integrand in the exponent 
at momentum scales below an infrared factorization scale k. We rewrite Eq. (7) as 
the sum of a perturbative term, labelled with the subscript PT, where all p± > k, 
and a soft term that contains all NP physics: 



In 



{u,Q,a) 



jf + Tf A (a , (pL)) J^(e---^r. l) 



dp] 



pi 



p±/Q 



du 



pI/Q 2 



In 



+ £>-term 

— cr PT (u, Q, k, a) 




^p n ± A(a s (p ± )). (8) 
PI 



We have suppressed terms of order 0(v/Q 



2-a 



~ a /Q 2 )i which include the entire 



S-term of Eq. (7), as indicated. Introducing the shape function as an expansion in 
powers of v/Q, from the expansion of the exponent in the second equality of (8), with 
NP coefficients A n (/c), we arrive at, 



(z/,Q,a) 



Op 



(v, Q, k, a) /o,np I ~q ,k 



(9) 
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1,1 Mi- K ) s T^| A "W("i) ■ (10) 

We find the simple result that the only dependence on a is through an overall factor 
1/(1 — a) which leads to the scaling rule for the shape function [10]: 

i 

l-a 

(11) 

For example, given the shape function for the thrust, a = 0, at a specific cm. energy 
Q, one can predict the shape function and thus from Eq. (9) the complete cross section 
including all leading power corrections for any other value of a. 





Figure 1: Differential distributions {1/ a tot )da /dr a for a = 0, and a = —0.5 at Q — 
91 GeV. Dash-dotted line: output of PYTHIA, a = 0; dotted line: perturbative 
contribution at NLL/NLO, a = 0, as defined in Eq. (8); dashed line: same as dotted 
line, at a = —0.5; solid line: prediction for a = —0.5 using Eqs. (9) and (11). The 
data are taken from Refs. [25, 26, 27]. 

This is illustrated in Fig. 1, where we plot the differential distributions for the 
thrust (a = 0) and for angularity a = —0.5. The dotted (a = 0) and dashed (a = 
—0.5) lines are the theoretical predictions at NLL from Eq. (7) matched to fixed 
next-to-leading order (NLO) calculations with EVENT2 [23]. In order to compute 
the shape function via Eq. (9) for a = 0, we use PYTHIA's [24] output at a = 
(dash-dotted line), instead of fitting a function to the thrust data, since PYTHIA fits 
these data very well. The prediction for the full differential distribution for a = —0.5 
(solid line) is then found via (9) and (11) from the so determined shape function and 
the NLL/NLO computation (dashed line). 
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The above derivation rests on two main assumptions: Our starting point, the 
NLL resummed cross section (7), describes independent radiation from the two pri- 
mary outgoing partons. Inter-hemisphere correlations are neglected, although they 
are present in the resummed formula valid to all logarithmic orders [8]. However, it 
has been found from numerical studies that such correlations may indeed be unim- 
portant [20, 21, 22]. Furthermore, we assume in the separation of perturbative and 
NP effects (8), that long-distance physics has the same properties under boosts as 
the short- distance radiation. Success or failure of the scaling (11) would thus provide 
information about these properties of long-distance dynamics. 

4 Non-Global Logarithms 

The parameter-dependence of the family of angularities can also help to control the 
sensitivity to wide-angle soft gluon radiation of certain less inclusive, but still infrared 
safe, so-called non-global, observables. 

4.1 Non-Global Observables 

Non-global observables measure the radiation into only part of phase space, an exam- 
ple is the energy flow into interjet regions [11, 28]. However, these observables retain 
sensitivity to radiation away from the location of interest via secondary radiation, as 
shown in Fig. 2 [29, 30, 31]. Fig. 2 illustrates an energy flow observable associated 
with radiation into a chosen interjet angular region, Q with the complement of Q 
denoted by Q. We are interested in the distribution of Qn for events with a fixed 
number of jets in Cl, 

A + B^ Jets +X n + R n (Qn). (12) 

Here Xq stands for radiation into the regions between Q and the jet axes, and Rq for 
radiation into Q. 

As shown in Fig. 2, for these kinds of observables there are two main sources 
of large logarithmic corrections: "primary" emissions, such as gluon 1 in Fig. 2, 
are emitted directly from the hard partons into Q. Phase space integrals for these 
emissions contribute single logarithms per loop, of order a" ln n (Q/<5n) = «™ln n 1/e, 
where we have introduced the variable e = Qq/Q. These logarithms exponentiate as 
above, (7), and may be resummed in a straightforward fashion [7]. There are also 
"secondary" emissions originating from the complementary region 0, illustrated by 
configuration 2 in Fig. 2. As emphasized in Refs. [29, 30], emissions into Q from 
such secondary partons can give logarithms of the form a™ {^(Qq/Qq), where Qq 
is the maximum energy of radiation in Q. These have become known as non-global 
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Figure 2: Sources of global and non-global logarithms in dijet events. Configuration 
1, a primary emission, is the source of global logarithms, configuration 2 results in 
non-global logarithms. 

logarithms. If no restriction is placed on the radiation into 0, then Qq can approach 
Q, and the non-global, secondary logarithms can become as important as the primary 
logarithms. The non-global logarithms arise because real and virtual enhancements 
associated with secondary emissions do not cancel each other fully at fixed Qn, or 
equivalently e. 

Such non-global logarithms cannot be factorized into a fixed number of jets as in 
Eq. (7), and therefore do not exponentiate straightforwardly. However, it was found 
in Ref. [32] that in the limit of large numbers of colors the radiative corrections to non- 
global energy flow observables can be described by a non-linear evolution equation, 
which is formally equal to the Kovchegov equation that describes the high-energy 
behavior of the ^-matrix [33]. Moreover, for a specific non-global observable, the 
heavy-quark multiplicity in a certain region of phase space, the non-global evolution 
linearizes [34, 35] and becomes formally equal to the BFKL equation [36, 37]. 

4.2 Event Shape/Energy Flow Correlations 

An alternative to resumming the secondary logarithms numerically or analytically 
in the limit or large numbers of colors, is to limit the radiation into the secondary 
part of phase space by studying correlations between the non-global logarithms and 
another global observable. Thus "mini-jets" of high energy Qq — > Q are suppressed, 
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because the radiation into is now restricted by the value of the global observable. 
If the global observable in addition has an adjustable parameter, such as the family 
of angularities, Eq. (5), the importance of the non-global components relative to the 
primary energy flow can be controlled. 

To be specific, we discuss the energy flow/event shape correlations for dijet events 
in e + e~ annihilation where color flow does not add extra complications. We will 
discuss color flow in hadronic collisions in the next section. We measure the energy 
flow into Q by weighing the final state iV with 

< N ) = k E<&> ( 13 ) 

where Qi denotes the energy of parton % which is emitted into the region of interest, 
Q. This is correlated with the event shape (5), measured in fl. Thus the cross section 
is given by 

M dedr a Q) = W V \ M W\ 2 8(e-e(N))6(r a -T a (N)) ( 14 ) 

where we sum over all final states N that contribute to the weighted event, and where 
M(N) denotes the corresponding amplitude for e + e~ — > N. 



4.2.1 Resummation 

We are again interested in dijet events, which correspond to the limit of small r a 
and e, r a ,e <C I. In this limit, as discussed above, the cross section receives large 
logarithmic corrections in both variables. However, since for the shape/flow correla- 
tion all of phase space is taken into account, we can apply standard factorization and 
resummation techniques [38, 39]. 

We arrive at the following factorized cross section in convolution form: 

' sQ e r as Q\ 



a tot de dr a 



S S 




5 (r a - r as - r aJl - T aj2 ) . (15) 



Here, \x is the factorization scale which we take for simplicity equal to the renomal- 
ization scale. H is a short-distance function, where all momenta are far off-shell, of 
order Q. This hard scattering function is therefore independent of the shape/flow 
correlations, and only depends on the energy (and implicitly on the scattering an- 
gles). S is a subdiagram that describes soft radiation away from the two jets Jj, 



8 



which contain all information about radiation collinear to the primary partons, and 
are thus constrained only by the event shape variable r a . In contrast, S depends on 
both, the event shape variable r tt , and the wide-angle radiation into Q measured via 
e. In Laplace transform space the convolution (15) becomes a simple product, 

doie^ = r dTae -, a r a da(e,r Q) 
de Jo de dr a 

1 da (e, v„ Q) = g ,Q\ g UQ _ ^ Q_\ JQ\ _ ^ 



Vtot de \fl "ai*/ iLip \ V aH 

where now the unbarred quantities denote the Laplace transforms of the corresponding 
functions in (15). 

Resummation of the large logarithms in e and u a is now straightforward. We use 
the fact that the physical cross section is independent of the factorization scale 

d ( da(e u a ,Q) \ = Q 
d/i \ de J 

to obtain the following renormalizat ion-group equations for the soft function, 

H-£-\nS ( — ,ev a ,— J = -isMfi)) , (19) 

an \ fl Vail) 

where 7 S is the soft anomalous dimension. Analogous equations are found for the jet 
functions. The anomalous dimensions can depend only on variables held in common 
between at least two of the functions. Because each function is infrared safe, while 
ultraviolet divergences are present only in virtual diagrams, the anomalous dimensions 
cannot depend on the parameters is a , e or a. This leaves as arguments of the 7s only 
the coupling a s {ji). The solution to Eq. (19) is found easily, 

JeQ Q\ Q (eQ Q \ ~ I **{a.W) 

S —,ei/ ai — =S —,ei/ 0) —,o e ™ . 20 

V A* V a.n J V A*0 Va^O J 

Choosing /i of order v a relegates all large logarithms lnz/ a into the exponent, as 
desired. By employing further evolution equations and proper choices for the ini- 
tial scales we achieve exponentiation of all large logarithms lnz/ a and lne, leaving 
only logarithms of order ln(eu a ) unexponentiated. For the full, a bit unwieldy final 
resummed expression and all technical details, we refer to Refs. [8, 40]. 

To summarize, large logarithms stemming from either the restriction of radia- 
tion by the global event shape or by the non-global observable can be resummed 
by standard techniques into a Sudakov form analogous to Eq. (7), leaving only an 
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unexponentiated piece from soft wide-angle radiation that contributes logarithms of 
order ln(eu a ), or equivalently, upon transformation back from moment space, of or- 
der ln(e/r a ). Moreover, we observe that the not yet fully exponentiated logarithms 
ln(e/r a ) stem from wide-angle soft radiation that completely decouples from the jets. 
This soft radiation can be described by functions constructed entirely out of non- 
abelian phases, or Wilson (eikonal) lines [38]. It is well-known that such eikonal 
quantities exponentiate directly by a reordering of color factors and with the help 
of an identity for eikonalized propagators [41, 42]. Therefore all large logarithms, 
including logarithms of {e/T a ) exponentiate, as was also found in Ref. [43] by an 
alternative analysis via the coherent branching formalism. Furthermore, by choos- 
ing a parameter-dependent global event shape we are able to dial the importance of 
these correlated logarithms ~ ln(e/r a ), and thus study the behavior of the non-global 
logarithms. 

5 Color Flow 

Above, we have discussed energy flow observables, which provide information that is 
in some sense complementary to what we learn from the study of event shapes. Event 
shapes shed light on the global distribution of radiation, with the main features de- 
termined by the underlying hard scattering. Long-distance effects soften these main 
features, but not drastically so. Interjet energy flow, on the other hand, reflects the 
interference between radiation from different jets [44], and encodes the mechanisms 
that neutralize color in the hadronization process. The study of the interplay be- 
tween energy and color flow in hadronic collisions [45] may therefore help identify the 
underlying event [46, 47], to distinguish QCD bremsstrahlung from signals of new 
physics. 

5.1 Hadronic Angularities 

In order to study color flow we need to define event shape observables suitable for 
hadronic collisions. Because parton-parton scattering is highly singular in the for- 
ward direction, with a 1/ sin 4 (#/2) behavior due to gluon exchange, e + e~ event shape 
definitions have to be adapted to accommodate this behavior in the hadronic case. 
For further discussion of hadronic event shapes we refer to Ref. [2] . Again, the intro- 
duction of a parameter-dependent observable allows us to control the aforementioned 
difficulties. 

We propose the following extension of the set of angularities to hadronic colli- 
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sions [48], 

r a (N) = -1= E P±,i e_l??<l(1 " a) . (21) 

Here, in contrast to Eq. (5), the transverse momentum p± ;i and rapidity r)i of parton 
% is measured relative to the collision axis in the cm. frame. The parameter a allows 
us to control the sensitivity to the forward direction. 

The hadronic cross section for the process A + B — > Jets + X at cm. energy y^s 
is now given by a convolution of standard parton distribution functions (PDFs) </>//# 
with a short-distance partonic cross section da: 



d@ ab(j~o,i P±) 
dr)dp±dr a 



E J dx A dx B <f>f A /A (xa, (J>f) <fif B /B (%b, Vf) 



xs( P± --^-) (22) 
\ 2 cosh 77 y dfjdr a 

li F denotes the factorization scale. The PDFs <f>f/ H describe the probability for finding 
parton / in hadron H with momentum fraction x. Hatted quantities are given in the 
partonic cm. frame, which can be found from the corresponding hadronic quantities 
via fj — r) — (1/2) \o.{xa/ 'xb), and s = x A x B s. The superscript (f) denotes the Born- 
level 2^2 processes, 

f: /U + /B-/1 + /2, (23) 

and the sum is over all possible processes. The hard scale is set by the transverse 
momentum of the observed jet, p±. Corrections to this leading-twist factorization 
begin in general with powers of h^ cu /p]_ due to multiple scatterings of partons. 

As above, the partonic cross section will receive large logarithmic corrections due 
to gluon radiation. However, in the case of (21), r a 7^ at lowest order, due to the con- 
tributions of the outgoing jets (recall that we measure with respect to the beam axis). 
Thus, we encounter large logarithms of order lnf a , f a = r a — J2j=i,2P±,J e ~^ r>J ^ 1 ~ a \ 
where the sum is over the outgoing jets. In what follows, moments are taken with 
respect to f a . 



5.2 Color Evolution 

The partonic cross section, da^ / (df)dr a ) can be refactorized analogous to Eq. (17), 
but now we have to take the non-trivial color flow due to colored initial-state partons 
into account. The regions that give leading contributions are as in the previous 
section a hard scattering, soft, and jet functions, two for the outgoing jets, and two 
jet functions for the beam jets. The latter jets must be defined to avoid double 
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counting due to the parton distribution functions. Furthermore, care must be taken 
in the definition of the outgoing jet functions due to the logarithmic enhancement 
proportional to lnf a instead of lnr a . Such definitions are quite non-trivial, and we do 
not attempt a full treatment here. Instead, we sketch the main features that emerge 
independently of the particular details of the chosen observable. 

As was observed, for example, in Refs. [11, 45, 49], there is no unique way of 
defining color exchange in a finite amount of time since gluons of any energy, including 
soft gluons, carry octet color charge. The functions from which we construct the 
refactorized partonic cross section are therefore described by matrices in the space of 
possible color exchanges. This is because as the (re)factorization scale changes, gluons 
that were included in the hard function become soft and vice versa, as illustrated in 
Fig. 3. Due to intrajet coherence [44], however, the evolution of the jets themselves is 
independent of the color exchanges. Once a jet is formed, collinear radiation cannot 
change its color structure. Therefore, the refactorized partonic shape/flow correlation 
can be written in moment space as 



da < - i \p a ,p ± ,ii F ) 
drj 



H { L ) I {p±,V,fi,fi F ) Sfl(u a ,fj,n) ]J Jc\va,P±,V,n,fi F ). (24) 



7 (f) 



analogous to the corresponding correlation in e + e~ events, Eq. (17). Now the hard 
and soft functions, H and S, respectively, are matrices in the space of color flow. 
Repeated indices in color space, L, I, are summed over. The superscripts (f ) label the 
underlying partonic process, as in (23). For a list of convenient color bases for the 
various 2 — > 2 processes see for example Refs. [40, 45]. fi is a refactorization scale, 
not necessarily equal to the factorization scale in Eq. (22). 




\i=>\i 




Figure 3: Cartoon illustrating color evolution. H and S are the hard scattering and 
the soft function, respectively, as described in the text. 



The color flow of the event shape is captured by the renormalization-group equa- 
tion (RGE) that results quite analogously to Eq. (19) from the requirement that the 
physical cross section be independent of the refactorization scale. Since H and S are 
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matrices in color space, the RGE is now a matrix equation, with anomalous dimension 
matrices T, 

-^a,^)(r%, a 3 (n)j) JL . (25) 

Since the anomalous dimension matrices for global event shapes are found as usual 
from virtual graphs only, the matrices are independent of the shape function and thus 
universal. They are tabulated in various places, for example in Refs. [40, 45], and 
have been implemented in the automated resummation program CAESAR [50]. 

For example, for the computation of the soft anomalous dimension matrix for the 
partonic subprocess qAQB — > <?i<?2 the t-channel singlet-octet basis is convenient, 



J_ 1 

2N A ' ri rB ' r2 2 



where r- label the color indices of parton i, and N c is the number of colors. The soft 
anomalous dimension matrix in this basis is given by, 

r (w-55) _ ( 2CpT ~% U \ (07) 

-\-2U -±-{T-2U) J' {27) 

where T and U are functions of the Mandelstam variables s, t, u in the partonic cm. 
frame: T = In [=fj +iir,U = In (^) + in. 

The partonic cross section resumming single logarithms in the color flow is then 
found by solving Eq. (25), 

da (f) (u a ,p±,fi F )) 



dfj 



II 4 i \va,P±,V,fi,fi F ) (28) 



A 7 

where we have transformed the matrices H and 5 to a basis, where the anomalous 
dimension matrices T are diagonal. 

(r2W)) 7/J - \fm„ = <> (r2W))„ R^il (29) 

with \y the eigenvalues of T, and R is the transformation matrix. Greek indices 
(3, 7 indicate that a matrix is evaluated in the basis where the eikonal anomalous 
dimension has been diagonalized. 
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From Eq. (28) we see that the color flow at NLL is encoded in event shape inde- 
pendent anomalous dimension matrices. It may therefore be possible to find a suitable 
range of the parameter a in Eq. (21) or perhaps another, different, event shape observ- 
able, where small-x effects and problems due to incomplete detector coverage in the 
forward direction become unimportant, and the study of color flow is unobstructed. 
Moreover, by combining the analysis of color flow presented above with the study of 
energy flow discussed in the previous section, as outlined in Ref . [40] , we may extract 
further important information about the dynamics of QCD. 

6 Conclusions and Outlook 

Above, we have illustrated how the introduction of a parameter-dependence can fur- 
ther increase the amount of information that can be extracted from the study of event 
shape observables. Aside from information about particle spin and various parame- 
ters of QCD, such as color factors and the running coupling, event shapes, especially 
their distributions, can be used to study the interplay of short- and long-distance 
dynamics, and obtain insight into the mechanisms of confinement, energy and color 
flow. 

The study of the family of angularities is but a starting point. Given the amount 
and versatility of event shape observables for processes ranging from electron-positron 
annihilation over deep inelastic scattering to hadronic collisions, there are many more 
possibilities and applications. It is certain that continued theoretical and experimen- 
tal efforts in this direction, at present and future accelerators, will reveal a wealth of 
information, and probably a few surprises. 
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